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Abstract 



Weak decays and mixing of S-hadrons play a special role in our understanding of 
the physics of the Standard Model and beyond. The measured amplitudes, however, re- 
sult from a complicated interplay of weak and strong interaction effects. Understanding 
strong interaction dynamics, which becomes simpler for heavy quarks, is thus a fun- 
damental part of the game. In this review, several theoretical aspects of B dynamics 
which are relevant for current and future experimental measurements are discussed. 



1 Introduction 



Hadronic weak decays are a unique source of information on the physics of the Standard 
Model and beyond. Experimental measurements and theoretical studies of these decays 
may provide us the key for understanding still unsolved problems such as the reason for the 
existence of several fermion families, the origin of fermion masses and weak couplings and 
the mechanisms which give rise to CP violation. Through quantum-loop effects, they are 
also a window which allows us to explore and constraint new physics effects at high energy 
scales still unaccessible to particle accelerators, such as those due to heavy supersymmetric 
particles. Unfortunately the physical amplitudes result from a complicated interplay of weak 
and strong interaction effects, which must be disentangled in order to extract the values of 
the fundamental parameters. 

In this framework, i?-hadron decays play a special role, since the 6-quark belongs to the 
heaviest, and less "coupled", of the three known families and its partner, the top quark, 
has a mass which is of the order of the weak scale, ~ Gp"^"^ . This observation, for 
example, has led to several speculations: that spontaneous symmetry breaking occurs via 
the Nambu-Jona Lasinio mechanism that quark masses are dynamically determined 
by radiative corrections or that the fundamental quark mass matrix, at large energy 
scales, is constrained to have very specific forms [Q]. Another aspect which is specific to 
5-hadron decays is that the mass of the 6-quark is much larger than the scale of strong 
interactions, m;, ^ ^qcd- In this situation, QCD dynamics becomes much simpler and in 
many cases, by using the Heavy Quark Effective Theory (HQET) [Q, we can obtain accurate 
theoretical predictions to compare with the experimental measurements. An example is given 
by the determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vcb\ from 
exclusive {B 0*1^1 and B —>■ Dli/i) and inclusive (h — > XJi^i) semi-leptonic decays. Since, 
however, the mass of the 6-quark is large, but not infinite, a related problem, to be discussed 
in the following, is that of controlling power corrections, i.e. corrections suppressed as inverse 
powers of the heavy quark mass. 

In this article, I review several theoretical issues in the physics of B decays which are rel- 
evant for current or future experimental measurements. A special attention is paid to the 
discussion of the systematic errors present in the different theoretical approaches, and of the 
possible progresses that we may expect in the near future. In preparing this talk, it has been 
unavoidable to make a selection of subjects. Among those for which new theoretical studies 
or new experimental measurements appeared this year, I have chosen those topics which 
look to me particularly interesting, and/or for which I have enough compentence to express 
a personal point of view. For lack of space, I am unable to cover here radiative decays and 
supersymmetric effects in FCNC processes. A discussion on these points can be found in the 
talks by A. Ali and A. Masiero at this workshop. For further information on other subjects 
related to 5-physics, not covered here, the reader can also refer to some recent review papers 
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by A. Ali §, A. Buras [|] and M. Neubert §. 

This article is organized as follows. In brief overview of leptonic decays, B^-B^ 

mixing and CP violation is given; in sec. 5-meson exclusive semi-leptonic decays, and the 
extraction of \Vch\ and \ Vuh\, are discussed; some aspects of exclusive non-leptonic decays are 
reviewed in sec. |; power corrections in the HQET are discussed in sec. || some considerations 
on inclusive semi-leptonic decays can be found in sec. |; the present status, and unsolved 
problems, for 5-hadron lifetimes will be summarized in sec. 0, followed by a brief conclusion, 
which contains also an outlook on possible improvements in the accuracy of the theoretical 
predictions, sec. ||. 



2 Z^, B^—B^ mixing and CP violation 



The leptonic decay width of the B meson is given by 

■^2 IT/ _ 12 / 

Mb 1 




A precise knowledge of the leptonic decay constant /s, analogous to /tt in vr — > decays Q, 
would allow an accurate extraction of the CKM matrix element \ Vub\i which is actually known 
with a relative error of about 25%, see below. One of the most precise determinations 
of fs is presently given by lattice QCD calculations and it is shown in table |I], together 
with other results for several physical quantities of interest for 5-physics |P]-|[T^. The 
error on fs reported in the table comes mainly from the extrapolation of the lattice results 
to the continuum limit and does not includes the systematic error due to the quenched 
approximation, which has been estimated to be of the order of 10-30% |Tl|. A global 



theoretical error of 25 MeV, corresponding to a relative error of about 14%, would provide 
us with an accurate determination of \Vub\i through a measurement totally independent from 
the semi-leptonic channel, and with experimental and theoretical systematic effects different 
from that case. Even though the rate (which goes as the squared lepton mass) is much 
smaller, the B^ /^^^'^ leptonic decay channel is also quite interesting, because the much 
smaller number of events could be compensated by a much higher efficiency in detecting the 
final state charged lepton. 

fs also enters in phenomenological analyses of the B^-B^ mixing amplitudes, together with 
the so called renormalization group invariant 5-parameter 13 b- The square of the mixing 

parameter ^ = fB\J Bb is in fact related to the matrix element of the renormalized Ai? = 2 
Hamiltonian [l^. All theoretical calculations of Bb tend to give values very close to one 



for both the B^ and the mesons, see table |I] and ref. |]T^ for QCD sum rules. Thus 
^ I use the normalization convention in which /^r — 132 MeV. 
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fo (MeV) 


/d. (MeV) 


Ib (MeV) 


Ibs/ Jb 


205 ± 15 


235 ± 15 


175 ±25 


1.15 ±0.05 


extrapolated 


Bb, hqet 


BbJBb, 


id (MeV) 


1.40 ±0.10 


1.08 ±0.06 ±0.08 


1.01 ±0.01 


207 ±30 



Table 1: Lattice predictions for several quantities of interest in heavy flavour physics. A 
detailed explanation of the methods used to compute these quantities, and of the evaluation 
of the systematic errors can he found in /|^-/|7^/. The numbers given here are my personal 
compilation of lattice results obtained from the numbers reported in refs. and the 

results which can be found in ref. / |7^/ . 



the strength of the mixing is essentially regulated by the meson decay constant. Taking 
into account the theoretical uncertainties, and the experimental errors on the relevant CKM 



matrix elements and on the B lifetime, one predicts ^Mb^ = 0.55 ± 0.13 ps |T5[, to be 
compared with the experimental world average AM^^ = 0.464 ± 0.018 ps. 

The value of is also important for CP violation in B decays. It was realized long time 
ago that a value of ^ ~ 200 MeV, combined with a large value of the top mass, (in 
the following I will use mf'^ {mf^^) = 167 ± 8 GeV) leads to a large value of sin 2/3, the 
parameter which controls CP violation in B ^ J/tpKg decays 

This is shown in fig. |I| where the distribution of acceptable values of sin 2/3, obtained from a 
combined analysis of the kaon CP violation parameter e and of the B^-B^ mixing rate AMb^ 
is shown [^. The histogram refers to the distribution of values of sin 2/5 obtained by varying 
within their errors the values of the different experimental and theoretical quantities which 
enter in the theoretical predictions. The dashed histogram is obtained as the previous one, by 
imposing the further condition that C, = 207 ±30 MeV. From the dashed distribution shown 
in the figure, one gets sin 2/5 = 0.68 ± 0.10, a value quite reassuring for future experiments 
at the 5-factories. 

The reliability of the lattice estimates of the meson decay constants is demonstrated by the 
fact that fos, which was predicted long before the first experimental measurement [T^, agrees 



with the present experimental world average from Dg /iz/^ decays, fn^ = 241 ± 21 ± 30 
MeV [fu^ ~ 260 MeV if one includes the L3 measurement from Dg tu^-) [|TB|. Further, 



indirect information can be extracted from the experimental lower bound on the Bg-B^ 
mixing amplitude. Taking AMb^ > 9.2 ps from experiments and fB^/fB = 1-15 ± 0.05"'"'^'^^ 
from theory, where the last error is an estimate of the error coming from the quenched 
approximation 0, one obtains > 160 MeV in good agreement with the direct theoretical 

^ In the ratio of the decay constants many systematic errors cancel. 
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Figure 1: Distributions of values for cos 6 and sm2j3, forrrit = (168±8) GeV. The values of 
the parameters are those used in the second of refs. /[7^. The solid histograms are obtained 
without using the information coming from the experimental determination of AMb^ ■ The 



dashed ones uses AMb^ = 0.464±0.018, with the condition 177 MeV<^d = fB^yBBa < 237 
MeV. 4 



prediction |TI 

I want to conclude this section with a comment on the "heavyness" of the charm quark. The 
HQET is not a model. It is QCD systematically expanded in powers of a small parameter 
fJ'QCo/^Q, where mg is the quark mass and fiqcD stands for a scale of the order of the 
strong interaction scale. Thus if /igcD = ^qcd ~ 300 MeV, since rricharm ~ 1.2-1.5 GeV, 
the expansion parameter is a small number even in the charm case. If, instead, figcD = 
2vr//) ~ 1.2 GeV or of the order of the nucleon mass, the expansion parameter is of order 
one and the heavy quark expansion is expected to fail. Lattice studies of the behaviour of 



the meson decay constant /m as a function of the inverse meson mass P]-[10[ have shown 



that the expansion parameter is, at least in this case, of the order of 1 GeV. These results 



have been confirmed by calculations using QCD sum rules [|14|. Such a large value of Hqcd 
implies that power corrections are ~ 80%-100% for the charmed meson decay constant, and 
~ 25%-30% in the i?-meson case. Corrections of the same order are also expected, and were 
indeed found in lattice calculations, for the semi-leptonic form factors relevant in B —>■ tt and 
B ^ p decays 0. A special, and probably unique, case where the application of the HQET 
to the charm is succesful, because of Luke's theorem [|l^, is represented by the semi-leptonic 
B ^ D and D* decays which are discussed in sec. ^. In all the other cases, as in the 
discussion of the charmed hadron lifetimes of sec. |^, there is strong experimental evidence 
that the charm is not heavy enough to make HQET applicable. 



3 Exclusive semi-leptonic decays 

In this section I will discuss some theoretical issues related to exclusive heavy-to-heavy and 
heavy-to-light semi-leptonic decays. 



3.1 Vcb from B D^vi 

In heavy-to- heavy decays, a golden plated determination of \Vcb\ is given by the comparison 
of the differential decay rate for the process B D*lvi to its theoretical expression ^ 

aw 48 TT'^ 

y^w + l {Ms-Mn*y j \^c,\:F[w), (3j 



For lack of space I will not discuss the process B Dlvi for which both measurements, and theoretical 
predictions, exist. 



5 



where the variable w, related to the momentum trasfer q, is given by 



W = Vr ■ Vn* = = ■ 4 

In the infinite mass limit, up to logarithmic corrections, J^{w) tends to the Isgur-Wise 
function ^{w). When both and mcharm become very large, the normalization of ^{w) at 
zero recoil is fixed by the conservation of the weak current, i.e. ^{w = 1) = 1. In these 
limits then, up to the trivial kinematical factors appearing in eq. (Q), a model independent 
extraction of \Vcb\ is possible. Two elements are important for an accurate determination of 
\Vcb\, however: the control of the corrections due to different effects at zero recoil, and the 
extrapolation of the rate to w = 1. 

i) The normalization of J-'{w = 1) is modified by short-distance strong interactions effects, 
which can in principle be computed in perturbation theory, and by non-perturbative l/mq 
corrections (starting at second order by virtue of Luke's theorem |]TP|). We have then 

J^{w = l) = r^Jl + C2^^ + ...)=r]A(l + 6^/^2) . (5) 

A recent two-loop calculation has reduced the uncertainty coming from the perturbative 
corrections to less than 1%, giving rjA = 0.960 ± 0.007 pOf. The evaluation of 5i/m2 is 
more problematic, since this is a non-perturbative quantity. Two different approaches were 
followed in the past. On the one hand, a classification of all the ^/mq local operators of 
the HQET contributing to Si/m^ was made, and their matrix elements estimated in a model 



dependent way These estimates gave —di/m^ = (3 ± 2)%. On the other hand, zero- 

recoil sum rules, which are in principle model- independent, gave a much larger correction 
corresponding to — 5i/m2 = 5.5%-8% |2^. After some controversy, M. Neubert, by using 



a "hybrid" method, found a "Solomonic" result, which is "not-in-disagreement" with both 
the values given above, namely —Si/m^ = (5.5 ± 2.5)% [^. This prediction is currently 
adopted in all the estimates of the theoretical error on \ Vcb\ extracted from the experimental 
measurements of the exclusive semi-leptonic decays. 

A recent study of the accuracy of power corrections in effective theories, illustrated by 
several explicit examples, showed that the predictions from the zero recoil sum rules, though 
model independent, are potentially subject to large corrections coming from higher order 
perturbative terms which have not been computed yet |2^. The same arguments apply to 



the calculations of refs. [p^, p3 . For this reason, and given the high accuracy reached by the 



experimental measurements, I believe that a further effort is needed to try to reduce, and to 
evaluate more realistically, the error on ^1/^2. This is very difficult because the deviations 
due to 5i/m2 are rather small, and no theoretical method, including lattice simulations, can 
control the uncertanties to a sufficient degree of accuracy. 
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ii) Since the rate vanishes atw = l (see eq. (|^)) in order to extract J^{w = 1) an extrapolation 
of the data is necessary. This extrapolation is usually done by expanding the form factor 
around w = 1 

J^(w) = ^(l)(l -p\w-l) + c{w -If + ...). (6) 

The slope p and the curvature c can be fitted directly from the experimental measurements 
of the differential rate |2^, It is, however, helpful for reducing the uncertainty of the 



extrapolation to have theoretical estimates of these parameters, for which several bounds, 
and calculations, obtained using different techniques, exist. The spread of the theoretical 
predictions for p is, unfortunately, still rather large while for c there are only analyticity 
bounds which relate its value to that of p pSf . In analogy to what is done experimentally, p 



has been extracted from lattice calculations by fitting the form factor computed in numerical 
simulations to eq. (P). This fitting procedure, however, introduces a rather large uncertainty 
in the final result . I want to mention that a method to obtain directly the p and c, without 
any fitting procedure, or perturbative calculation, has been proposed in ref. but only a 
preliminary feasibility study of this technique has been done to date |3^. Contrary to the 
case of Si/m'2, it is likely that more accurate theoretical results for p will be obtained in the 
near future. 

As shown in fig. |^, the values of J^{w = l)|Kb| extracted by the different experiments are 
quite consistent. This is true in spite of the fact that the fitted values of p are somehow 
different From the world average of J-'{w = l)|V^fc|, by assuming from the theory 

J^(w = 1) = 0.91 ± 0.03 one obtains, 

iKfel = 0.0382 ±0.0018 ±0.0012 , (7) 

where the last error comes from the estimated theoretical uncertainty. 



3.2 Vub from B plui and B nlvi 



For heavy-to-light semi-leptonic decays, HQET is less useful since the normalization of the 
current is not fixed by symmetry. Still scaling laws can be established for 
example, up to logarithmic corrections, one finds pT| , 



AMy^ 



1 ± 



5 



+ ... 



(8) 



where fsiPn) is the form factor relevant in B ^ n decays. P denotes the heavy pseudoscalar 
meson, Mp its mass, the pion momentum and A and 6 are non perturbative quantities. 

^ Following ref. [Q, I use here p to denote the physical slope, in order to distinguish it from p which refers 
to the slope of the Isgur-Wise function ^. 
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Figure 2: Determinations of \Vcb\^{w — 1) from different experimental measurements. The 
average is also given. 
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The scaling law in eq. and others which can be derived for the form factors relevant in 
B — >■ plui OT B ^ K*'y decays, are valid in the heavy meson rest frame, for small values 
of the light meson three-momentum, i.e. for i4 Mp. Scaling laws at ~ can also be 



derived by computing the form factors using light-cone wave functions . For example one 
finds Ai{q'^ = 0) ~ Mp^^"^, where Ai is the axial form factor entering in i? ^ plui decays. 
Unfortunately, corrections to the asymptotic scaling behaviour were found to be large both 



in lattice calculations ||^, ^ and using QCD sum rules [Q, cf. also the discussion on power 
corrections in sec. @. A further constraint which can be used in theoretical predictions is 
that coming from the equality of some of the form factors in the helicity basis at = 0, e.g. 
f^{q^ = 0) = f^{q'^ = 0), and from analyticity The above information is particularly 
useful for lattice calculations, where only a modest range of the physical q^ is presently 
accessible. 

Several theoretical methods have been used to predict the form factors for B ^ n and B ^ p 



decays. Besides the lattice calculations and the QCD sum rules [§3| mentioned above, 
several quark model calculations have also appeared in the literature. A complete list of 
theoretical calculations can be found in refs. [1^, |2^, |2^. It is interesting that, with the 
increasing experimental precision, it is possible to discriminate the "good" from the "bad" 
models by confronting the ratio of widths r{B —>■ p)/T{B — >■ vr), in which the unknown \Vub\ 
cancels, to the theoretical predictions ^ A compilation of results, obtained by the 
CLEO collaboration using different theoretical models is shown in fig. ^. In the extraction 
of \Vub\i some models, the results of which are given in the figure above the horizontal line, 
have also been used in the Monte Carlo analysis of the data, in order to generate events to be 
confronted to the experimental measurements. The corresponding values for \Vub\ are only 
given in these cases. The last value, denoted as "Endpoint method" refers, for comparison, 
to the result obtained from inclusive charmless decays, to be discussed in the following. From 
the figure, we see that the dominant error is the theoretical one. This is true in spite of the 
encouraging progresses made by the theoretical calculations in the last few years. These 
progresses are evident if one compares the results of fig. ^ to those available at the time 
of the 5th Symposium on Heavy Flavour Physics, held in Montreal in 1993 |^. A further 
improvement of the present situation is likely to occur in the near future, particularly from 
lattice calculations. 



4 Exclusive non-leptonic decays 

In this section, I briefly discuss a few items connected to exclusive two-body non leptonic 
decays. 

Factorization : Factorization is an interesting phenomenological idea which allows us to pre- 
dict non-leptonic amplitudes in terms of semi-leptonic form factors and of vector and pseu- 
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Figure 3: Experimental determinations of \Vub\ from exclusive S — > tt and B ^ p semi- 

leptonic decays. The results depend on the theoretical model used in the experimental analysis. 
Only the models above the orizontal line have been used to extract \Vub\ from the data. The 
value of \Vub\ from inclusive decays is also given at the bottom of the figure. 
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doscalar decay constants. In a way, it corresponds to an elaborated version of the vacuum 
saturation approximation. Factorization was originally used to estimate kaon decay ampli- 
tudes In the eighties, it was applied by Stech and collaborators in analyzing D-meson 
decays |^ and its explanation in terms of the expansion, where N is the number of 



colours, was subsequently proposed |3^. More recently, it has been suggested that fac- 



torization holds in the heavy quark limit, for some energetic decay channels, because the 



ultrarelativistic light mesons in the final state do not have time to exchange gluons [p9, 40 



In ref. |^T| it was shown, however, that the demonstration of factorization given in has 



some problems. It may be still true that factorization becomes valid in the infinite mass limit, 
but a convincing proof of this intuitive physical idea is still missing. Several observations 
are necessary at this point. 

i) Factorization does not work for charm meson decays, where, moreover, large final-state 
interaction effects are measured (final-state interactions are zero if amplitudes factorize). 
This may be due either to a failure of factorization, or to the fact that, as observed above, 
the charm quark is not heavy enough. Indeed factorization, combined with the assumption 
that final-state interactions are dominated by nearby resonances, has been succesfully applied 
to phenomenological studies of charmed meson decays, for example in ref. . 



ii) The effective Hamiltonian which controls weak charm decays has the form 

^ = EC.(^)o,(/i), (9) 

where Oj(/i) are composite operators renormalized at the scale /i and Cj {Mw / jj) the corre- 
sponding Wilson coefficients. The dependence on /i of the operators is compensated by the 
dependence on /i of the Wilson coefficients. The factorized matrix elements of the operators, 
however, are //-independent. Thus the final physical amplitude, in the factorization approx- 
imation, seems to depend on the renormalization scale. It may be true that ther exists a 
scale at which non-factorizable contributions are small, e.g. fi = nib for 5-decays, but a 
violation of factorization must be present to some extent. 



iii) As recently shown in ref. [^] factorization and quark-hadron duality are incompatible. 
To reconcile duality and factorization a certain amount of non-factorizable contributions has 
to be allowed. 

In conclusion, we know that non-factorizable contributions ought to be there, and that they 
are important for charm decays. The basic question is then to have a quantitative estimate 
of their size in i?-decays. This problem has been studied in extensive experimental analyses 
of two-body exclusive decays and it could be possibly accessible to lattice simulations 
in the near future see below. 



Non-leptonic decays on the lattice: There exist only quite old (and not very accurate) studies 
of heavy-quark two-body non-leptonic decays on the lattice, performed at ttiq ~ mcharm 
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DE 



CE 




Figure 4: Non-penguin Feynman diagrams for non-leptonic two-body decays. 



15[. In these calculations only "non-penguin" amplitudes were considered, namely those 
corresponding to the following Feynman diagrams: disconnected emission (DE), denoted 
also as T or T' in ref. the colour-suppressed connected (non-factorizable) emission (CE), 



denoted also as C or C; the hehcity suppressed disconnected annihilation (DA), denoted 
also as A or A', and connected annihilation (CA), denoted also as E or E' . These diagrams 
are given in fig. ^. 

In it was found that the amplitude of the non-factorizable diagram CE was about one 
third of the leading contribution given by DE, as expected on the basis the expansion, 
and that the amplitude of DA was indeed compatible with zero. The amplitude correspond- 
ing to CA, which is assumed to be negligible in some theoretical analyses [Q, was found, 
however, to be as large as one half of that from DE. All the above results were unfortunately 
afflicted by large statistical and systematical errors, and have never been repeated since. 
The reason for this is that all lattice calculations of two-body non-leptonic amplitudes were 
abruptly interrupted when the Maiani- Testa no-go theorem was published |^^. This theo- 



rem states the following. In numerical simulations on the lattice, it is unavoidable to make 
a Wick rotation to the Euclidean space-time, where in- and out-states are indistinguishable. 
From the standard study of the correlation functions one obtains then the average of the 
in- and out- amplitudes instead of the physical one, thus loosing all information about final 
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state interactions Q. For example, in the calculation of the decay of a 5-meson into two light 
mesons Mi and M2, mediated by the weak Hamiltonian 7^^, one gets 

= ^( out{M, M2\n'^\B) (Ml M.in'^lB)) , (10) 

instead of out{Mi M2\TC^\B) only. Moreover, at large time- distances, the correlation func- 
tions from which the amplitudes are usually extracted are dominated by the off-shell matrix 
element (Mi(pi = 0) M2(p2 = 0)\7{^\B), instead of the physical amplitude. In general, 
the Maiani- Testa theorem prevents the calculation of time-like electro-magnetic form factors 
and strong and weak decays in two or more hadron final states. Its extension can be used 
to show that it also limits the possibility of computing inclusive cross-sections and decay 
rates [H9|. A possible way out of this problem, under mild assumptions, although at the 



price of introducing some model dependence, has been recently suggested in ref. We 



do not know at present if the proposal of ref. MM will, in practice, succeed in extracting the 



physical amplitudes including the phases due to final state interactions. Feasibility studies 
are currently underway. If the proposal works, a new field of investigation (factorization, 
scaling laws, effects due to nearby resonances) will be open to lattice calculations. 

Other interesting two-body amplitudes: 

i) Penguin Pollution: in B^ tt+vt" decays, CP violating effects are given by the combined 



contribution of tree-level diagrams and penguin diagrams |£0|- order to extract the rel- 
evant phases, it is the necessary to discriminate the two effects. Experimentally, this can 
be achieved by separating the term proportional to cosAM^^t from that proportional to 
sin AMs^t in the time-dependent asymmetry 

_ iV(i?° ^ 7r+7r-)(t) - iV(^g ^ 7r+7r-)(t) 

iV(i?o^7r+7r-)(t)+iV(50^7r+7r-)(t)- ^ > 

A theoretical calculation of the tree-level and penguin contributions using the method of 
ref. P3, combined with the experimental analysis, would certainly help to reduce the sys- 



0^0 



tematic uncertanties in the final results. 

ii) Isospin analysis: by combining the measured rates of the B^ tt+tt^ and B^ 
decays, one can remove the penguin contribution, by isolating the pure A J = 3/2 ampli- 
tude |5^. The difficulty of this method is due to the measurement of the decay rate into 
two neutral pions, which is expected to be colour suppressed. An alternative method, based 
on a combined study of B^ tttt, 5" — ^ ttK and B^ — > KK decays, using approximate 
flavour SU{3) symmetry has also been proposed. There is no space to discuss the details. 



and uncertainties, of this method, for which the reader can refer to I^Ol . I want only to add 



one observation, which may stimulate further investigation. In refs. 0], it is usually 
^In ref. p9], a similar problem was found in the context of QCD sum rules. 
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assumed that the non-f actor izable annihilation diagram CA is neghgible, in contrast 
within large errors, and for smaller heavy quark masses, with the results of ref. [45 



but 



see 



above. A reappraisal of the lattice calculations of the non-leptonic amplitudes, on large lat- 
tices, and by using the theoretical improvements developed in the last few years |^ would 
certainly be very useful in this respect. 



5 Power corrections in the HQET 

In this section the problem of controlling power corrections in the heavy quark effective theory 
(HQET) is discussed. These corrections are necessary for an accurate evaluation of both 
inclusive and exclusive decays of hadrons containing a 6-quark. The basic question is whether 
they can be computed with a theoretical error smaller than the corrections themselves |25 



To be specific, we start by considering the generic expression for the inclusive decay rate 



m-m 



6A2 „5 , 167r^/iM. 



and the ratio of form factors in B ^ D*lvi decays 0] 

h^{uj = vv') C^{uo) / A ^, , \ 

i = rxH ^ + — + • • • • ^3 

Ai^(a; = V ■ v') C^[uj) \ mc J 

In the above expressions are the relevant vector and axial vector form factors; Cj, Cj, 
Cf, C^iuj) and C^(u;) are short distance quantities that can be computed in perturbation 
theory; the quark mass mi,, the binding energy A, the kinetic energy Ai, the chromo- magnetic 
term A2 and the -Bf 's, that are related to the matrix elements of four-fermion operators, 
require instead a non-perturbative evaluation. 

In some cases the non-perturbative parameters can be directly related, at the leading or- 
der in l/mq, to a physical quantity. This is the case of the chromo-magnetic term which 
corresponds to the B*-B mass splitting, A2 = 1/4(M^* — M^). In general, however, the 
coefficients of the l/mq expansion, and the corresponding non-perturbative parameters, are 
plagued by renormalon ambiguities. This problem is present for example in the definition 
the binding energy A appearing in eq. (^) and in the matrix elements Bf of the four fermion 



operators in eq. (|T2]) [|^, [^-[^. In these cases, the cancellation of the ambiguities occurs 
between the matrix elements of the higher dimensional operators and the Wilson coefficient 
functions (C|, Cj, ...). The latter have to be computed to sufficiently high order of perturba- 
tion theory for the cancellation to be under control. Let us imagine that we want to achieve 
an accuracy of 0(l/m^). A simple argument indicates that, in principle, the leading Wilson 
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coefficient has to be computed at least to an order such that as{mb)'^ ~ {AQco/^b)'' ||2^, |K 
Since, in the calculations performed up to now, only the ffist few terms of the perturbation 
series are known, it is not possible to check that the remaining terms are negligible. By 
studying some simple examples, it has been shown that the knowledge of only a few terms 
is, in general, insufficient to control the power corrections . Although the results rely on 



some approximations that were adopted in these examples, it is likely that the conclusions 
will remain valid in general, and that in most further theoretical effort is needed. 

Let us consider, for example, the unphysical parameter A = Mh — fripoie which is frequently 
used in phenomenological studies of 5-physics. Results and bounds for A have been presented 
in the literature |^-|^. In order for this to make any sense, A must be defined precisely 
in terms of some physical quantity, and its value will depend on this physical quantity and 
on the order of perturbation theory used to extract A. In ref. it has been shown that, 
even by using a consistent and renormalon-free definition of the binding energy, one expects 
that the value of A will change by several hundred MeV, i.e. by an amount comparable to 
its value, as the order of perturbation theory is increased. The error on A can be seen as an 
error on the MS mass of the 6-quark m^^(m^^) as derived from the HQET ||6^ 



m'^^m'^'') = 4.17 ±0.05 ± 0.20 GeV , (14) 

where the last (and larger) error comes from an estimate of the uncertainty due to higher 
order corrections in perturbation theory. 

Similar problems are encountered in the case of Ai. One example is given by the zero recoil 
sum rules [^, In this case, in order to derive bounds on the relevant form factors, 

and on the heavy-quark kinetic energy parameter Ai, a physical cut-off A is introduced, to 
suppress the contribution of states with excitation energies greater than A. The difficulty 
in achieving an accurate determination or bound with a physical cut-off is due to the higher 
order perturbative corrections and only by arriving to an order such that a" < 0(Ai/A^) 
it will be possible to obtain significant results. A connected example is provided by the 
different theoretical determinations of Ai, which are given in table Some of them are 
not compatible, if we accept as realistic the errors quoted in the original papers. I believe 
that the differences are due to higher order perturbative corrections, which have not been 
computed, or estimated, yet. 



6 Semi-leptonic inclusive decays 

The idea that inclusive decay rates of heavy quarks can be computed in the parton model 
is quite old and was used, for example, in ref. to predict charm hadron lifetimes. To 
account for bound state effects, the partonic calculation was successively improved by the 
introduction of a phenomenological model, called the "spectator model" |7^, |75[, which has 
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Theoretical estimates of Ai 



Reference 



Method 



-Ai (GeV2 



Eletsky and Shuryak | 
Ball and Braun ^ 

Bigi et al. [0 
Kapustin et al. 



Ligeti and Nir pU| 



Gimenez et al. |70|] 
Chernyak |]7T | 



Gremm et al. F 



Neubert |7| 



QCD sum rules 
QCD sum rules 
ZR sum rules 
ZR sum rules +0{as) + 0{a^j3o 
Experiment 
Experiment 

Lattice 
Experiment 
Experiment 
Virial Theorem 



0.18 ±0.16 
0.52 ±0.12 
> 0.36 

< 0.63 if A > 240 MeV 

< 0.10 if A > 500 MeV 

-0.09 ±0.14 
0.14 ±0.03 
0.19 ±0.10 
0.10 ±0.05 



Table 2: Some of the values of Ai obtained in different theoretical analyses. "Experiment" 
denotes the extraction of Xi from the experimental data, for example the charged-lepton 
spectrum distribution in semi-leptonic B -meson decays. "ZR sum rules" denotes the zero 
recoil sum rules. 



been and continues to be extensively used to extract \Vch\ and \Vuh\ from inclusive semi- 
leptonic decays. A noble theoretical framework for the spectator model was then provided 
by the operator product expansion applied to inclusive decays of heavy quarks ||5l|l-||55[|. 
The general result is summarized in eq. ([T^), which holds for both inclusive semi-leptonic 
and non-leptonic decays. The first term on the r.h.s. of eq. (p!2D is the parton model result; 
corrections of 0(l/mf) are due to the kinetic energy — Ai and the spin interaction A2 of the 
heavy quark in the hadron 0; the four-fermion operator contributions to the 1 / correction, 
written in terms of the Bf parameters, correspond, in the old language, to W-exchange and 
interference higher-twist effects. These effects were considered in order to explain the lifetime 
differences among the different charmed hadrons ||76|| . 

In modern language, the use of the parton model to compute inclusive decays relies on the 
assumption of "quark-hadron" duality: in semi-leptonic decays, where the integration over 
the lepton and neutrino momenta corresponds to an integration over the invariant mass 
of the recoiling hadronic system, we have global duality, whereas in non-leptonic decays, 
where the hadronic mass is fixed to be equal to the heavy hadron mass, we speak of local 
duality. Duality is assumed, but cannot be proven in QCD. We know that it provides a good 
description of several physical processes such as inclusive deep inelastic scattering, inclusive 
e"'"e~ annihilation into hadrons, Drell-Yan processes, etc. The assumption of global duality 
seems also to work for semi-leptonic (and radiative) decays. This means that we are able to 

® Up to a redefinition of tlie quarlc mass, there is indeed a relation between tlie Fermi momentum of the 
heavy quark in the model of ref. [ [zs) and Ai. 
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predict the inclusive TsL^b XJi'i), Tsiib XJi^i) and TsL^b s'j) rates. This is not 
the case for local duality since there is strong experimental evidence that we are unable to 
explain 5-hadron lifetimes in the framework of the l/m^ expansion, see sec. ^ 

The situation for TsL^b XJui) is substantially unchanged since last year. For this reason, 
I direct the reader to ref. for a discussion of the theoretical uncertainties in the calculation 
of this quantity. Here I limit myself giving the upgraded world average values of \Vcb\ from 
inclusive decays 

= 0.0391 ± 0.0007 ± 0.0030 , 
\Vcb\^^^ = 0.0412 ± 0.0004 ± 0.0030 , (15) 

where the results from low energy and LEP energies have been given separated. The reason 
is that the different central values are induced mostly by the differences in the semi-leptonic 
branching ratio measured in the two cases ||18|. Notice the excellent agreement between 
these results and the number given in eq. (|^. 

The above results make us confident that we are able to make accurate predictions for 
^SLib XJui). A related problem is whether we are also able to predict the semi-leptonic 
branching ratio Bsl, and the average number of charmed hadrons ric- Since the former is 
given by the ratio of TsL{b XJiyi) to the total width, the problem is whether we are able 
to predict inclusive non-leptonic widths. Recently, there was a lot of debate about whether 
the theoretical value of Bsl is in agreement, or in disagreement, with the experimental 



measurements [p5| , [77| , [78| and, in the second case, several explanations, within and beyond 
the standard model, have been proposed [^-[Q. Since there are important differences 
between low energy and LEP results 



= 10.49 ± 0.17 ± 0.43 , = 1.12 ±0.05 

= 10.95 ± 0.13 ± 0.29 , n™ = 1.22 ±0.08 (16) 

it is probably premature to draw any conclusion at this stage. 

\Vub\ could be extracted, in principle, from the total semi-leptonic width in the same way 
as \Vcb\- In practice, however, only the endpoint spectrum is used, in order to eliminate 
the background due to b ^ dui decays. The cut on the lepton energy Ei is fixed by the 
requirement that the recoiling hadronic system X has an invariant mass Mx < Md, so as 
to inhibit charm decays. In this region, however, we run into troubles with the operator 
product expansion. Schematically, we can write 



dT{b — > ului) 
dEi 




Tniy) + ■■■ , (17) 



where . . . represent terms which are suppressed as powers of rrib. Only the first few terms in 
the expansion above are known because they can be related to Ai, A2, etc. The higher order 
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terms depend on matrix elements of higher dimensional local operators, which can at most 
be estimated within a specific phenomenological model. 

Near the end point, 

ymin ~ 1 — ^x/^t {.Umax IS always of 0(1)) SO that when Mx ~ 



AqcdMb, all the terms in the expansion ( |T7| ) become equally important. At even larger 
values of Ei, Mx ~ ^qcd- In this case we are in the elastic region (corresponding to B ^ p, 
-B — s> vr, etc.) and the decay cannot be analyzed with operator product expansion. Unfor- 



tunately, Mx ~ y AqcdMb is about 1.2-1.3 GeV, a value very close to the kinematical cut 
which is currently used in the experiments to extract \Vub\- This implies that the extraction 
of \Vub\ is still subject to an intrinsic model dependence, which can be reduced either by 
enlarging the region of Ei from where \Vub\ is extracted, at the price of increasing the charm 
background, or by a sound theoretical prediction of the "wave" function ||5^, Q of the heavy 
quark. 



7 Non-leptonic inclusive decays 

In this section it will be shown that it is very unlikely that higher order corrections in l/rrib 
can explain the discrepancy between the theoretical predictions for the B to Ab lifetime ratio 

R = ^b/ta^- 

For R, the dependence on the fifth power of the quark mass, as well as all renormalon 
ambiguities up to order I /ml, cancel out. We thus expect R to be one of the soundest 
predictions of HQET. The discrepancy between theory and experimental data is evidence 
for the failure of local duality in the calculation of the non-leptonic (NL) widths of heavy 
hadrons. 

To order 1/ml in the heavy quark expansion, the lifetime ratio R is given by 

2m( ml 

where the difference of kinetic energies \i{B) — Ai(At) = 0.01 ± 0.03 GeV^ can be obtained 
from the spectrum of the heavy hadrons; A2(A;,) = 0; \2{B) = (M|* - M|)/4 = 0.12 GeV^ 
and Cg/3 ~ 1.1. Using the numerical values given above one obtains R = 1.02 + 0(l/mf), 
incompatible with the observed ratio of lifetimes R = 1.27 ± 0.05. Within the standard 
HQET, in order to explain the large deviation from unity for this ratio it is necessary the 
the matrix elements of the four-quark operators entering in the l/mf corrections (and in 
particular those relative to the A^) be much larger than those estimated in quark models, 
and more recently with QCD sum rules |80|. In my opinion, this remains an open, but 
very unlikely possibility. An alternative explanation is that for NL widths there is evidence, 
contrary to the orthodox expectations, of corrections of 0{l/mb). These corrections are well 
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described by the simple ansatz that replaces the quark mass with the decaying hadron mass 
in the mq factor in front of the NL width [§, 7^ 



r^2 ™5 i\/r5 



This replacement provides a much better description of the NL widths. In ref. ||7^, it was 
shown that, for beauty, both the problems of the SL branching ratio and of the difference 
in the lifetimes of the baryon and the B mesons are quantitatively solved. For charm 
a much better fit to the seven known lifetimes is obtained in terms of four parameters of 
reasonable size: one lifetime, one interference contribution for D^, one for S+ and a smaller 
W-exchange term for D'^. 1 will now give more details on the analysis of ref. |]7D| . 

Given the agreement between the value of \Vcb\ extracted from inclusive b ^ c and exclusive 
B D* semi-leptonic (SL) decays and the approximate equality of the SL widths for 
charmed hadrons (within large uncertainties and with a rather limited number of cases 
where measurements exist), it is assumed that the HQET is able to predict the SL widths 
TsL of the heavy hadrons. Under this assumption, and using the ansatz (|19|), one gets 

R = [1 _ 2.24BsLm + 2.2ABsl{B) + 0{l/Ml) , (20) 

where the factor 2.24 arises from taking the electron, the muon and the tau SL modes in the 
ratio 1:1:0.24 [Q, and the difference in the SL rates of the A;, and B has been neglected. From 



Ma, = 5623 ± 6 MeV [|1 and Mb = 5279 ± 2 MeV [||], by using Bsl{B) = 10.77 ± 0.43%, 
one finds f\ 

— = 1.29 ±0.05 (21) 
in perfect agreement with the experimental result. 

So far this agreement can be considered to be a numerical accident, lacking theoretical 
explanation. For this reason, it is quite interesting to consider what happens in the case 
of charmed particles for which seven lifetimes have been measured, see table |^. At lowest 
order in the l/mg expansion, a much better agreement with the experimental results for the 
lifetimes is obtained by replacing the heavy-quark mass by the hadron masses as in eq. ([191) . 
We neglect at this stage any other mass correction and write Tj^l^Mh) = TtotiMn) — '^^^sl, 
where for a universal value chosen as the average of the experimental values for D^, 
and Ac, or Tsl = (0.174 ± 0.015)ps^^ is used The dependence on the hadron mass of 
Tnl{Mh) will be taken according to Tnl{Mh) = {Mh / Mq)'^T n 1(^0) with n = 5, where Mq 
is around the average mass of the relevant hadron. We then have 

TUMh) = r-\Mn) = ^"'(^o)(^)" + 2Tsl{1 - (^)") (22) 
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Hadron 


Mass (MeV/c2) 


r (ps) 


Bsl{%) 


TsL = Bsl/t (ps-^) 




1869.4 ±0.4 


1.057 ±0.015 


17.2 ± 1.9 


16.3 ± 1.8 




1864.6 ±0.5 


0.415 ±0.004 


8.1 ± 1.1 


19.5 ±2.6 




1968.5 ±0.7 


0.467 ±0.017 






K 


2285.1 ±0.6 


0.200 ±0.011 


4.5 ± 1.7 


22.5 ±8.5 


' — 'c 


2470.3 ± 1.8 


0.098 ±0.019 






'=:± 


2465.1 ± 1.6 


0.350 ±0.055 








2704 ± 4 


0.055 ±0.023 







Table 3: Properties of charmed mesons and baryons; the VL^ values are an average of the 
data quoted in ref. j^^. 
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0.8 — 




M (GeV/c^) 



Figure 5: Lifetime vs. mass for charmed particles. The dashed line is the best fit described in 
the text for all seven points assuming proportionality of the NL widths to Mfj, Mh being the 
hadron mass. The solid line is the best fit restricted to only the D^, A^, H° and Qc lifetimes 
with the same assumptions as before. 
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Let us fix n = 5, Mq = 2.3 GeV and Tsl = 0.174 ps~^ and fit all seven known lifetimes in 
terms of r(Mo). T(mo) = 0.181 ps is then obtained. The corresponding fit is shown in fig. ^ 
(dashed curve). We see that four out of seven lifetimes are in very good agreement with the 
fitted curve. The lifetimes of D^, of H"*" and, to a lesser extent, of Dg are clearly out. The 
discrepancies for and S+ can be attributed to the interference effect Note that 
is the only meson that can have interference at the Cabibbo allowed level and S"*" is the only 
baryon that can have double interference, in the sense that S"*" = cus and both u and s can 
interfere with the corresponding quarks from c sud. For Dg the observed smaller difference 
is attributed to the possibility of W-exchange fS^. All of these effects are of order f^/ml or 



1/m^. The solid line in fig. |^ has been obtained from a modified fit where only the D^, A^, S*^ 
and flc lifetimes have been considered. In this case, t(Mo) = 0.161 ps, with the respectable 
value of the /d.o.f. given by ~ 3.5, is obtained. For comparison, the fit of the quark mass 
to constant lifetimes results in a jd.o.j. ~ 250. Finally, for the same four lifetimes, one 
fits the power n in eq. (p^ , keeping fixed the value of Mq and r(Mo) at the observed values 
for the meson. In this way one may check whether the best power for n is close to 5. 
n = 4.5 ± 0.5 is found, where the error arises from the experimental errors on the lifetimes. 



1 - (-^] 



where 



Moreover, if we write for S"*" and Ds the expression r ^ = r ^{Mh) 
T is the experimental number given in table |^ and t{Mh) is taken from the previous fit to 
the four remaining lifetimes (with n=5 ), /i = 1.6, 2.2 and 1.3 GeV are obtained for D^, S+ 
and Ds, respectively. We see that the resulting values of this correction are large, as it is 
obvious from fig. |^, but not unreasonable. 

The conclusion is that a number of experimental facts show that Tnl for charm and beauty 
decay approximately scale with the fifth power of hadron masses (apart from corrections of 
order l/rriQ or smaller) rather than having ©(l/m^) corrections to the universal behaviour 
predicted by the HQET. Since M^^ = (mt + Aji/J^ = m^(l + 5Aj:/^/mb) + . . ., this implies that 
the corrections are of 0{l/mf,). Since the validity of the operator expansion in the timelike 
region, in the vicinity of the physical cut, is not at all guaranteed ||88|, the failure of the 
short distance approach can be attributed to a violation of the local duality property that has 
to be assumed for NL widths. This is to be compared to the case of the SL widths for which 
the predictions of the HQET are not inconsistent with the data. The experimental evidence 
for NL widths calls for a reexamination of the underlying theoretical framework. An attempt 
in this direction can be found in ref . [^] . More extended and accurate measurements of the 
SL widths of each hadron separately are also needed in order to assess the validity of the 
heavy quark expansion in semi-leptonic decays. 



^ The same numbers as in ref. M| are used here since the analysis has not been repeated since. 
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8 Conclusion and outlook 



Beauty hadrons provide an extraordinary laboratory to test the Standard Model, to look 
for hint of new physics and to understand QCD dynamics. For this reason an exceptional 
experimental and theoretical effort has been invested in it. In this review, I have tried to point 
out those dynamical aspects which we ought to understand and control at a quantitative 
level in order to extract the relevant information on some fundamental parameters, such as 
the mass of the b quark, the relevant CKM matrix elements and CP violation. 

Much theoretical progress, and improvement in understanding, has occured this year in the 
framework of the HQET and of the calculation of inclusive rates. Still a lot of work remains 
to be done, particularly for those quantities where the theoretical uncertainties are larger 
than the experimental errors. Progress can be expected quite soon for the slope and the 
curvature of the form factors in i? — > 0*1^1 decays, in the understanding of factorization 
and in the calculation of non-leptonic two-body decay amplitudes. It appears more difficult 
to reduce the error in the normalization of the form factor at zero recoil for B D*lvi and 
to eliminate the model dependence in the extraction of \Vub\ from inclusive semi-leptonic 
decays. A real puzzle, which demands a deep reappraisal of current ideas, and prejudices, 
remains that of the A;, lifetime, which seems to indicate a failure of local duality in inclusive 
non-leptonic decays. 
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